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ABSTRACT 

Stereo data collected by the HiRes experiment over a six year period are examined for large- 
scale anisotropy related to the inhomogeneous distribution of matter in the nearby Universe. We 
consider the generic case of small cosmic-ray deflections and a large number of sources tracing the 
matter distribution. In this matter tracer model the expected cosmic ray flux depends essentially 
on a single free parameter, the typical deflection angle 9 S . We find that the HiRes data with 
threshold energies of 40 EeV and 57 EeV are incompatible with the matter tracer model at a 
95% confidence level unless 9 S > 10° and are compatible with an isotropic flux. The data set 
above 10 EeV is compatible with both the matter tracer model and an isotropic flux. 
1. Introduction 
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The observation of the cutoff in the spectrum 
of Ultra-High Energ y Cosmic Rays (UHEC Rs) 



!()()8al _ 

predic ted bv lGreisenl ( 1966t ): IZatsepin and Kuzminl 



(jAbbasi et alJ l2008at lAbraham et al. 2008 aj) as 



(|1966l ) provides compelling evidence for the short- 
ening of the UHECR propagation length at high 
energies. The highest energy events then must 
have come from relatively close sources (within 
250 Mpc). At these length scales the matter in 
the Universe is distributed inhomogeneously, be- 
ing organized into clusters and superclusters. One 
should, therefore, expect the flux of highest-energy 
cosmic rays to be anisotropic. 

In astrophysical scenarios, it is natural to as- 
sume that the number of sources within 250 Mpc 
is large, and that these sources trace the distri- 



1 



bution of matter. Under these assumptions, the 
anisotropy at Earth depends only on the nature 
and size of UHECR deflections. Measurement of 
the anisotropy, therefore, gives direct experimen- 
tal access to parameters that determine the deflec- 
tions, notably to the UHECR charge composition 
and cosmic magnetic fields. 

Several investigations of anisotropy in ar- 
rival directions of UHECRs have been previ- 
ously undertaken. At small angular scales, 
correlations with different classes of putative 



sources were claimed (e.g. Gorbunov et al.l 2004; 



Abbasi et all l2006t lAbraham et al.l 120071 l2008bl ) 



At larger angular scales and energies below 10 
EeV possible anisotro py towards the G al actic c en- 



ter was reported i n Havashida et al. ( 1999al lb) 



Bellido et al.l (1200 ll). but no t supported by more 
recent studie s (ISantosI 2007 ). At higher energies, 



Stanev et al.1 (| 19951 ) found evidence against an 
isotropic flux above 40 EeV through correlations 
with the supergalactic plan e, but this was not con- 
firmed by other authors (Havashida et al. 1996; 
Kewlev et aDll996t bird et al.lll998l ). Finally, us- 
ing the Pierre Auge r Obse rvatory (PAO) data, 
Kashti and Waxmanl ( 20081) have found correla- 
tions between UHECR arrival directions and the 
large-scale structure of the Universe which are in- 
compatible with an isotrop ic flux (see, however, 



compati ble witn an isotrop 
Koers and Tinvakovl[2009a| ) . 



In this paper, we analyze the data accumulated 
by the HiRes experiment for anisotropy associ- 
ated with the large-scale structure of the Universe. 
The H iRes experiment has been described previ- 
ously (lAbu-Zavvad et aDll999t iBover et all 12002 : 
Hanlonl 12008V It studied ultrahigh energy cos- 
mic rays from 10 172 eV to 10 20 2 eV using the 
fluorescence technique. HiRes operated two fluo- 
rescence detectors located atop desert mountains 
in west-central Utah. The data set used in this 
study consists of events observed by both detec- 
tors, analyzed jointly in what is commonly called 
"stereo mode" . In this mode the angular resolu- 
tion in cosmic rays' pointing directions is about 
0.8 degrees, and the energy resolution is about 
10%. The HiRes experiment operated in stereo 
mode between December, 1999, and April, 2006. 
At the highest energies HiRes has the largest data 
set in the Northern hemisphere. Large number of 
events, good angular resolution (better than the 
bending angles expected from Galactic and ex- 



tragalactic magnetic fields) and the wide energy 
range covered make the HiRes data particularly 
suitable for anisotropy studies. The exact data 
set used in this study was described previously in 
Abbasi et all (|2008bf ). 



We consider here a generic model that assumes 
many sources within 250 Mpc tracing the distri- 
bution of matter, which we refer to as the "matter 
tracer" model. We also assume that deflections 
of UHECR do not exceed the angular size of the 
nearby structures, that is 10-20°. In this regime, 
both regular and random deflections in magnetic 
fields can be modeled with a one-parameter dis- 
tribution, for which we take a Gaussian distribu- 
tion centered at zero with width 6 S . This width 
is treated as a free parameter, whose value we 
aim to constrain from the data. Constraints on 
9 S may then be used to obtain information on the 
strength of Galactic and extragalactic magnetic 
fields. In keeping with our assumption of small de- 
flections, we assume a proton composition in this 
study, which is consistent with the X max analysis 
based on the same dataset (for confirmation see 



I Abbasi et aTll2009h . 



The HiRes data is compared to model predic- 
tions using the "flux samplin g" test put forward by 



Koers and Tinvakovl (|2009at ). This test has good 



discrimination power at small statistics and is in- 
sensitive to the details of deflections. The compar- 
ison is performed at three different threshold en- 
ergies that have been used in previous studies: 10 
EeV, 40 EeV and 5 7 EeV (lHavashida et al.lll996l 



Abbasi et al.l 120061 : lAbraham et al.l 120071 ). An 



priori significance of 5%, corresponding to a con- 
fidence level (CL) of 95%, is chosen for this work. 

The paper is organized as follows. In section [5] 
we discuss the modeling of UHECR arrival direc- 
tions. Section |3] concerns the HiRes data used 
in the analysis, while section [4] describes the flux 
sampling method. We present our results in sec- 
tion [5] and conclude in section [5] 

2. Modeling of UHECR arrival directions 

Galaxy catalog — The distribution of matter in 
the local Universe is modeled with the 2 Micron 
All-Sky Redshift Survey (2MRS; iHuchra et al" 



2009) galaxy catalog, using galaxies as samplers of 
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the underlying matter density field0 The 2MRS 
is a flux-limited sample of galaxies, that is, the 
sample containing all galaxies with observed mag- 
nitude in the K s band m < 11.25. It contains 
spcctroscopically measured redshifts for all but a 
few galaxies. A number of cuts have been applied 
to the galaxy sample. First, the Galactic plane, 
where the sample is incomplete, has been excluded 
from the sample by removing all galaxies with 
|6| < 10°. Second, objects with D < 5 Mpc are 
removed because such objects should be treated 
on an individual basis0 Finally, the catalog is 
cut at 250 Mpc because the sample becomes too 
sparse. The resulting sample provides an accurate 
statistical description at smearing angles 9 S > 2°. 
The flux from sources beyond 250 Mpc is taken 
to be isotropic. A total of 15508 galaxies remain 
in the HiRes field of view after the cuts. To com- 
pensate for observational selection effects in the 
(flux-limited) 2MRS catalog, weights uf at are as- 
signed to the g alaxies with the sliding-box m ethod 
as described in lKoers and Tinvakov (2009b). 

Energy loss — UHECR fluxes are affected by 
energy loss due to redshift and interactions with 
the Cosmic Microwave Background (CMB). To ac- 
count for the resulting flux suppression, the inte- 
gral flux, ipi, from a single source is expressed as 
follows: 



<Pi(E,Di) 



J(E)f(E,D l 
4ttL>„ 2 



(1) 



where E is the threshold energy, D t is the source 
distance, J stands for the integral injection spec- 
trum at the source, and / represents the flux frac- 
tion that remains after interactions and redshift. 
We take an injection spectrum J(E) oc E~ 12 ex- 
tending to very high energies. The function / is 
determined using a numerical propagation code 
based on the continuous loss approximation tha t 
is described in lKoers and Tinvakovl (|2008l . l2009al ) . 
Energy loss due to interactions with the extra- 
galactic background light is neglected. In FigurcQ] 
top panel, the fraction / is shown as a function 
of distance for the different energies considered in 
this work. 



iThis sample was kindly provided by John Huchra. 

2 This corresponds to the ad hoc assumption that there are 

no UHECR sources within 5 Mpc. Different analyses are 

more appropriate to test this possibility. 



The strength of the isotropic flux component 
that is added to account for sources beyond 250 
Mpc also depends on UHECR energy loss. Us- 
ing the computer code described in the previous 
paragraph, we estimate the fraction g of total flux 
contributed by sources within 250 Mpc to be 0.4, 
0.7, and 1.0 for threshold energies E = 10 EeV, 
40 EeV, and 57 EeV, respectively (see Figure 1, 
bottom panel). 

Deflections — UHECR protons (as well as nu- 
clei) arc deflected by Galactic and intergalactic 
magnetic fields. These deflections are taken into 
account by an angular smearing procedure, which 
replaces the point-source flux, ip, by a flux density 
distribution: 



(fi ->■ ifi w s (9i) . 



(2) 



where w B (9i) represents the probability density 
that an UHECR is deflected by 9i, the angle be- 
tween galaxy i and the line of sight. This proce- 
dure also accounts for the detector's angular res- 
olution and prevents unphysical fluctuations due 
to the use of a catalog of point sources. In the 
absence of detailed knowledge on the structure 
of Galactic and extragalactic magnetic fields, we 
adopt a simple Gaussian probability density dis- 
tribution with characteristic angle, 9 S - This angle 
is treated as a free model parameter. The Gaus- 
sian distribution is a fair approximation when the 
deflections are small. For large deflections, details 
on the structure of the Galactic and extragalactic 
magnetic fields become important. Accounting for 
these details goes beyond the scope of the present 
study. 

Exposure — The HiRes exposure is modeled 
using our Monte Carlo simulation of the exper- 
iment (|Abbasi et all l2007t iBergmanl l2007t ) . The 
aim of this simulation was to create a set of Monte 
Carlo events that would be in all essences iden- 
tical to the actual data. In making the sim- 
ulation we put in the properties of cosmic ray 
air showers as measured by previous exper i ments 



(Bird et al 



993 lAbu-Zavvad et all l2000l l200ll 
lAbbasi et al. 20051 ). We used cosmic ray showers 
generated by the Corsika and QGSJet programs 



([Heck et al.lll998t iKalmvkov et alJll997l ) and sim- 
ulated the g eneration of fluores cence light (see 



references in Abbasi et al. 2008al) and its propa 



gation through the atmosphere (see references in 
Abbasi et al.l2008ah . A complete simulation of the 
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Fig. 1. — Top panel: Fraction / of integral CR flux 
that survives after interactions with the CMB and 
cosmological redshift function of distance D 
for threshold energies 10 EeV (solid line), 40 EeV 
(dashed), and 57 EeV (dotted). Bottom panel: 
Fraction g of total flux that is produced by sources 
within 250 Mpc, as a function of energy E. 



optics and electronics (trigger and front-end elec- 
tronics) of our detectors was performed. The re- 
sult was an excellent simulation of our experiment 
as evidenced by the very good agreement between 
data and simulated events in the distribution of 
all kinematic variables, e.g. zenith angle, impact 
parameter to detector, etc. By assigning Monte 
Carlo events times of occurrence taken from the 
actual on-time of the experiment we are able to 
calculate the exposure on the sky very accurately. 

Model flux maps — The integral UHECR flux 
from a given direction is expressed as follows: 



(3) 



where i enumerates galaxies in the 2MRS sample, 
w1 &t denotes the weight assigned to galaxy i in the 
catalog, w s (9i) is the deflection probability distri- 
bution, and <f>i SO is the UHECR flux arising from 
sources beyond 250 Mpc. 

The probability to observe a CR from a given 
direction is proportional to the product of flux and 
exposure. We denote this probability as 



$ = $3. 



(4) 



where S stands^for exposure. In Figure [2j the 
distribution of $ over the sky is shown for three 
different threshold energies. The contrast in the 
flux distributions becomes more pronounced with 
increasing energy. Also shown are the arrival di- 
rections of UHECRs in the HiRes data to which 
the model flux has to be compared. 

3. Data 

The data set used in th i s study was described 
previously in lAbbasi et all (|2008bt ). including se- 
lection criteria and a correction to the energy 
scale. Our sample of the 2MRS catalog does 
not cover the region near the Galactic plane with 
\b\ < 10°. We therefore removed cosmic ray events 
with |6| < 10° from the analysis. The resulting 
sample contains: 

• 309 events with E > 10 EeV; 

• 27 events with E > 40 EeV; 

• 10 events with E > 57 EeV. 

The arrival directions of these events are shown as 
white dots in Fig. [2] 
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4. Statistical test 







Fig. 2. — Hammer projection (galactic coordi- 
nates) of $ (flux times exposure) with threshold 
energies 10 EeV (top panel), 40 EeV (middle), and 
57 EeV (bottom). Darker gray indicates a higher 
value; the bands are chosen such that each band 
contains 1/5 of the total flux (weighted with ex- 
posure). Excluded regions, viz. the galactic plane 
(\b\ < 10°) and the region outside the HiRes field 
of view, are shown in white. White dots indicate 
HiRes events. All maps are produced with 6$ = 6°. 



The compatibility of a model flux map with 
the set of UHECR arrival directions is quanti- 
fied by the flux sampling m ethod introduced by 
Koers and Tinvakovl (|2009ah . The idea of the 
method is as follows. To any set of arrival di- 
rections one associates the set of flux values that 
are obtained by sampling a given flux map (such 
as the map shown in fig. [5]), i.e. by extracting the 
flux values at corresponding points on the sphere. 
The two-dimensional distribution of arrival direc- 
tions thus translates into a one-dimensional dis- 
tribution of flux values. If the reference model is 
true, this flux distribution will tend to high values 
since events fall preferentially into regions where 
the model flux (times exposure) is high. If, on the 
other hand, the reference model is not true, the 
flux distribution is more uniform because the cor- 
relation between arrival directions and regions of 
high model flux is (partly) destroyed. By compar- 
ing the flux distribution to a model flux distribu- 
tion, the compatibility between a set of arrival di- 
rections and model predictions can be quantified. 
This comparison is performed by the Kolmogorov- 
Smirnov (KS) test, which yields a test statistic 
D. The relevant statistical quantities, in partic- 
ular powers and p-values, are computed from the 
distribution of this test statistic. Note that this 
test does not involve any additional parameters 
like bin size. 

The ability of the test to discriminate between 
models may be quantified in terms of the statisti- 
cal power, P, i.e. the probability to rule out, at a 
given confidence level, the reference model when 
an alternative model is true. Within numerical 
uncertainties, the statistical power is equal to the 
fraction of event sets generated under the alter- 
native model that leads to rejecting the reference 
model. Figure |3] shows the number of events re- 
quired for a power P = 0.5 (i.e., a 50% proba- 
bility) to rule out (at 95% CL) the matter tracer 
model when the true flux is isotropic. The num- 
ber of events increases with increasing smearing 
angle and decreasing energy: the decreasing flux 
contrasts in the matter tracer model call for an 
increase in statistics to achieve the same discrim- 
inatory power. Observe that the event numbers 
indicated in Figure [3] are of the same order as the 
data analyzed in this work. We thus expect that 
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there is sufficient data to obtain meaningfui con- 
straints at 95% CL. 

5. Results 

5.1. Scan over smearing angles 

The ievei of compatibility between data and 
model predictions is quantified by a p- value, which 
represents the model probability of obtaining a 
measurement that is at least as extreme as the 
actual measurement. With our a priori choice of 
significance, a p- value smaller than 0.05 rules out 
the model. The probability of falsely excluding 
the model is then 5%, translating to a CL of 95%. 

Figure |4] shows the p- values obtained by the flux 
sampling method for the HiRes data and predic- 
tions of the matter tracer model. The smearing 
angle, 9 S , is treated as a free parameter. That is, 
at each value of 9 S and each threshold energy a 
flux map is generated and compared to the HiRes 
data as described above. The results can be sum- 
marized as follows: 

(a) For the threshold energies of 40 EeV and 
57 EeV, the tests show disagreement be- 
tween data and the matter tracer model for 
#s < 10° • Within this parameter range, a 
source distribution tracing the distribution 
of matter is excluded at a 95% confidence 
level. 

(b) For the threshold energy of 10 EeV, the test 
shows agreement between data and the mat- 
ter tracer model. 

The incompatibility between data and matter 
tracer model is illustrated by the non-correlation 
between the observed arrival directions and re- 
gions of high model flux shown in the two lower 
panels of Figure [5J 

We have also tested the data for compatibility 
with an isotropic model flux and found no dis- 
agreement, at 95% CL, for any of the three tested 
threshold energies (the data with threshold energy 
57 EeV are marginally consistent with an isotropic 
flux). 

5.2. Case study: E = 57 EeV, 9 S = 3.2° 

At energy threshold E > 57 EeV a corre- 
lation between the arrival directions of UHE- 
CRs and the location of AGNs contained in the 



12 edition of the Veron- Cett y & Veron catalog 



([Veron-Cettv and Veron 



the PAO ( Abraham et al 



20061) was r eported by 
20071 l2008bh . This cor- 
relation was found to be maximal for ip = 3.2°, 
where tji denotes the maximum angular distance 
between UHECRs and AGNs. In the Northern 
hemisphere, correlation with AGN was not con- 



firmed by the HiRes experiment (jAbbasi et al 
2008bh . 



Since AGNs are tracers of the distribution of 
matter in the Universe, the PAO result is sugges- 
tive of a more general correlation between UHE- 
CRs and the local structure of the Universe on an 
angular scale of a few degrees. The methods pre- 
sented in this paper allow a check on the existence 
of such correlations in the HiRes data. 

The results presented in the previous section 
disfavor a correlation between UHECRs and the 
local structure of the Universe. In fact, the flux 
sampling test yields p- values smaller than 10 -2 for 
the matter tracer model with 9 S < 6°, with a p- 
value of 7 x 10" 4 for 9 S = 3.2°. (Note that 9 S is 
not in 1 : 1 correspondence with if); both quantities 
are however representative of the angular scale of 
the problem). Focusing on the case of 9 S = 3.2° 
in more detail, Figure [5] shows the distribution of 
the test statistic D for the matter tracer model 
and for an isotropic flux for this smearing angle. 
The vertical line shows the value D bs = 0.59 ob- 
tained for the HiRes data. This demonstrates the 
strong incompatiblity between HiRes data and the 
matter tracer model for smearing angle 9 S = 3.2° 
and threshold energy E = 57 EeV. 

6. Conclusions 

To summarize, we have confronted the stereo 
data collected by the HiRes experiment with pre- 
dictions of the matter tracer model, a generic 
model of cosmic ray origin and propagation. The 
model assumes a large number of cosmic ray 
sources within 250 Mpc whose distribution traces 
that of matter, and relatively small deflections 
characterized by a single parameter, the typical 
deflection angle 9 S . We have found that the HiRes 
data with energy thresholds E = 40 EeV and 
E = 57 EeV are incompatible with the matter 
tracer model for 9 S < 10° at 95% CL. With an 
energy threshold E = 10 EeV the HiRes data are 
compatible with the matter tracer model. At all 
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three energy thresholds, the data are compatible 
with an isotropic flux at 95% CL. 

In the present analysis we have treated the de- 
flections as random and Gaussian, which is only 
appropriate for small deflection angles and lim- 
ited number of events. The actual deflections are 
expected to contain a coherent component due to 
the Galactic magnetic field. With the accumula- 
tion of UHECR events by PAO in the Southern 
hemisphere and by Telescope Array in the North- 
ern hemisphere, our analysis will become sensi- 
tive to the nature of deflections and thus, with 
proper modifications of the statistical procedure, 
may give direct access to the parameters of cosmic 
magnetic fields. 
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Fig. 3. — Number of events required for a 50% 
probability to rule out, at 95% CL, the matter 
tracer model if the true flux is isotropic. 
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Fig. 4. — The dependence of p-value indicating 
the level of (in) compatibility between HiRes data 
and model predictions on the smearing angle 9 S . 
Solid lines indicate a p-value equal to 0.05, be- 
low which the model is ruled out at 95% CL. The 
points represent numerical results (with estimated 
uncertainties of 20%); the lines are smooth inter- 
polations between these points. Top panel: data 
vs. matter-tracer model; bottom panel: data vs. 
isotropic distribution. 
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Fig. 5. — Model distribution of test statistic D 
when testing the matter tracer model, for both 
the matter tracer model ("Structure") and an 
isotropic flux distribution ( "Isotropy" ) . Here E = 
57 EeV and 9 S = 3.2°. The vertical line indicates 
the observed value £> bs = 0.59. 
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